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ABSTRACT 

The baro tropic vortlcity equation is utilized to investigate 
the steady topographic disturbance of an initially uniform zonal 
current. Solutions are obtained on the beta-plane for various dis- 
continuous bathymetric profiles. It is shown that a smooth slope 
does not appreciably modify the results. Stability criteria are 
derived that define the values of the flow parameters for which 
the streamlines become discontinuous. 

The solutions are evaluated for several observed disturbances 
of oceanic currents. The velocities implied by the observed wave- 
lengths are physically realistic, but the discrepancy between 
observed and calculated amplitudes suggests that the vertical 
structure of the velocity profile and the baroclinicity of the 
ocean are important. 
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INTRODUCTION 

In most studies of the large scale circulation patterns of the 
oceans, the sea floor is assumed to be either flat or infinitely deep. 

Such simplifications are made to facilitate the investigation of certain 
classes of flows in which it is suspected that topographic effects are 
nonexistent or of minor importance. That currents in the real ocean are 
to some extent topographically affected nas been pointed out by many 
authors. Perhaps the most universally cited example of such an effect is 
the wave-like structure in the dynamic topography of the Antarctic 
Circumpolar Current (ACC) where it crosses ►each of several ridges (see 
Figure 1 ) . 

Many contributions to the understanding of oceanic currents have been 
achieved through the application of the techniques and principles of 
dynamic meteorology. That many apparent similarities exist between oceanic 
and atmospheric flows was graphically demonstrated in 1951 when Fuglister 
and Wortliington published the results of an intensive multiple ship survey 
of a Gulf Stream meander. Their illustration of the motion and incipient 
occlusion of the meander bears a striking resemblance to an earlier 
illustration (Berggren, et al., 19U9) of the unstable waves observed at 
the 500 mb level in tne atmosphere. It will therefore be found instructive 
to consider briefly the results of previous investigations of topographic 
effects on both oceanic and atmospheric flows. 
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I - 1 . Oceanographic Investigations 

Ekman (1923), in an attempt to account for the deflection of the 
Gulf Stream over the Southeast Newfoundland Ridge, considered the effect 

A. * 

of bathymetric relief on a homogeneous current under a balance of pressure 
gradient, bottom friction, and constant Coriolis forces. His results show 
the current deflection to be dependent on the slope of the bottom, but not 
the total depth, and give a net lateral displacement of the streamlines (to 
the right for flow over a ridge in the Northern Hemisphere ) . For flow over 
successive ridges, Ekman's streamlines are 90 degrees out of phase with the 
topography, the maximum deflection occurring at the point of inflection of 
the bottom profile. 

In 1932 Ekman reconsidered the problem and included accelerations 
in addition to bottom friction, constant Coriolis force, and pressure 
gradients. Subject to the requirements that both the bottom slope and the 
variation in depth be small, he obtained a solution that gives streamlines 
in phase with the bathymetry, the maximum displacement to the right occurring 
over the point of minimum depth. 

In an attempt to rigorously investigate the effect of bottom friction 
on a topographically disturbed homogeneous current, GOrtler (19Ul) utilized 
Prandtl's boundary layer theory. He simplified the problem by considering 
a constant Coriolis parameter and a variation in depth that is small 
compared to the total depth of the fluid. His solution to the linearized 
equations reduces to that of Ekman (1932) provided that Gdrtler's friction 
parameters, including the boundary layer thickness, are properly 
identified with Ekman's "depth of frictional resistance." Gttrtler extended 
Ekman's solution by including the possibility of the upstream flow 
approaching the topographic feature at other than normal incidence. He also 
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evaluated the frictional effects for several physically reasonable values 
of the flow parameters and thereby demonstrated that such effects are 
negligible for most atmospheric and oceanic currents. Gflrtler pointed out 
though, as did Ekman (1932), that topographic effects in which inertia 
plays an important role may not properly be treated as independent of 
planetary effects when the scale of motion is such that the latter mast be 
considered. This point is central to the investigation at hand. 

Sverdrup (19Ul) took a new approach to the problem of topographic 
effects. He neglected friction and considered that the effect of topographic 
features was to disturb the internal distribution of mass. Being unable to 
treat the problem rigorously, he nevertheless was able to demonstrate 
qualitatively that a bathymetric feature which rises from a layer of no 
motion into an internal solenoidal field and its associated current will 
disturb that current in such a way that the physical sea surface will be 
depressed above the feature. This effect causes the current to deviate to 
the right (in the Northern Hemisphere), the maximum deviation being over 
the point of minimum depth, apd then to return to its original configuration 

I 

downstream. That the converse is true for flow over a depression is not 
clear, as there is no mechanism for disturbing the solenoidal field. 

Sverdrup's conclusions agree qualitatively with the observed features of 
the circulation in the Southern Ocean. 

Shtokman (19U7), in a discussion of his subsequently published (19ii8b) 
theory of topographic effects on ocean currents, demonstrated an apparent 
serious discrepancy between the theories of Ekman, GOrtler, and Sverdrup 
and observed flow features of the South Atlantic and Caspian Sea. Considering 
Defant's (19U1) chart of the "absolute" topography of the surface in the 
Atlantic, and relating this chart to the bathymetry of the South Atlantic 
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between 25° and 50° South, Shtokman observed that the eastward currents 
are deflected poleward over shoaling water and equatorward over deepening 
water. This is contrary to the observed deflection of the ACC. Since 
the direction of deflection in the earlier theories depended directly on 
the sign of the Coriolis parameter, there was no mechanism available to 
explain opposite deflections in the same hemisphere . 

Noting the above apparent discrepancy between theory and observation, 
Shtokman inferred that the Coriolis parameter, if considered constant, 
should not be the primary consideration. On this basis, he adapted an 
earlier (l9U8a) derivation of the equations relating wind stress to lateral 
friction to the problem of a variable depth. For zero wind stress, the 
signs of the terms in hi6 equations give qualitative agreement with 
Sverdrup's (19U1) conclusions. Shtokman's approximation of the complete 
equations including wind stress implies that the direction of the current 
deflection over a topographic feature depends on the sign of the local 
curl of the wind stress and not on the hemisphere. Thus, for shoaling 
water, a negative wind stress curl yields a deflection to the right, and 
vice versa. 

As Shtokman pointed out, his results agree with the observations he 
cited. Thus the absence of perturbations on Defant's 19Ul topography of 
the surface of the South Atlantic between IiO° and 50° South is explained 
by the existence of maximum wind stress, and hence zero wind stress curl, 
in these latitudes. A later presentation by Defant (1961) of the dynamic 
topography of the same area indicates substantially the same structure 
except for the addition of an anticyclonic eddy at approximately li5°S - 
20°E (see Figure 2). If this is indeed an area of zero wind stress curl, 
then Shtokman's theory does not account for this eddy. 
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The writer believes that a more fundamental objection to Shtokman's 
theory may be made on the basis of his method of introduction of a variable 
depth. The depth which Shtokman ultimately considers as representing the 
actual bathymetry was originally defined as that depth at which the 
transport components vanish. Charney (1955) and Fofonoff (1962) have 
shown that this depth of no motion is a function of the transport itself, 
and hence is of a nature that is fundamentally different from the 
bathymetry. 

Shtokman's refutation of the theories of Ekman (1923), Gbrtler (19U1), 
and Sverdrup (19U1) based upon consideration of the near shore currents in 
the Caspian Sea seems inappropriate in that those earlier theories specifi- 
cally excluded changes of depth that amount to a significant fraction of 
the total depth. The earlier theories also excluded flows initially 
parallel to the isobaths, such as are considered by Shtokman (19^7, Fig. 5)« 

Neumann (I960) attempted to combine analytically the results of 
Ekman (1923) and Sverdrup (19Ul). Neumann initially considered the steady 
flow in a baroclinic ocean to be maintained by a balance of Coriolis, 
friction (both wind stress and "internal" friction), and pressure gradient 
forces. He pointed out, as have Defant (l9iil) and others, that foi 
homogeneous water, the variation of Coriolis parameter with latitude requires 
that all purely geostrophic currents be zonal, unless the depth increases 
from the equator toward the poles in proportion to the sine of the latitude. 
Neumann obtained approximate solutions for several classes of flow, including 
that of homogeneous water flowing over a ridge under the influence of a 
variable Coriolis parameter. This latter solution gives a Sverdrup- type 
deflection of the streamlines without the consideration of a solenoidal 
field. Neumann's solutions are based upon the requirement that the change 
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of depth be small compared to the total depth, bat the numerical examples 
given do not clearly meet that criterion. 

Saint-Guily (1962) has extended his earlier (1959) treatment of the 
generalized Ekman problem to the case of a variable depth. Subject to the 
assumptions of vanishing frictioh and moderate changes in depth, Saint- 
Guily obtained solutions for the steady zonal flow of homogeneous water 
over a sill that correspond to those of Neumann (i960) for a variable 
Coriolis parameter. Saint-Guily' s initial disregard of the field 
acceleration terms in the equations of motion restricts the validity of 
his solutions to regions where the ratio of the relative vorticity to the 
Coriolis parameter (Rossby Number) is small. 

1-2. Meteorological Investigations 

As was noted above, topographic effects on ocean currents have 
parallels in meteorological problems. Queney (19U9) summarized previous 
theoretical treatments of topographic effects, including the results of 
nis own extensive treatment of the problem (I9li7) utilizing perturbation 
techniques. Queney considered the stability associated with the density 
stratification of an infinitely deep atmosphere, and it is therefore 
difficult to apply his results directly to a homogeneous ocean. Charney 
and Eliassen (19li9) concluded that the real gravitationally stable atmosphere 
may be treated as an "equivalent barotropic atmosphere" whose properties 
coincide with those of the real atmosphere at a height of zero horizontal 
divergence (approximately 500-700 mb). Charney and Eliassen were able to 
substantially improve their numerical prediction techniques by including 
a rough model of the effects of continental topography. The importance of 
such continental effects in the study of large scale atmospheric flows has 
been discussed by Bolin (1950). He noted the regular perturbations of the 
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pressure surfaces that occur on monthly (and longer) mean weather charts of 
both the Northern and Southern Hemispheres and argued in favor of a topo- 
graphic mechanism for the generation of such perturbations as opposed to a 
thermodynamic mechanism. Bolin linearized the equations and obtained 
approximate solutions for flow over an isolated mountain. 

Oi (1956) considered the problem of the frictionless barotropic flow 
of air over a mountain range under the influence of a constant Coriolis 
parameter. He was able to include accelerations in a manner similar to 
Ekman (1932) and, again for small changes in the total depth of the fluid 
column, obtained approximate solutions showing the effect of a mountain 
range on an initially zonal flow. 

1-3 . Initial Evaluation of the Problem 

All of the above investigations, with the exception of that of 
Shtokman, have essentially linearized the problem by considering only small 
changes in relief relative to the total depth of the fluid column. This 
may be shown (Bolin, 1950) to be equivalent to assuming that the relative 
vorticity of the fluid is small compared to the local value of the Coriolis 
parameter. In the ocean, with a mean depth of about li km, and in the 
equivalent homogeneous atmosphere, with a height of 8 km, it is apparent 
that there are many areas where such a linearization vail be inappropriate . 
At low latitudes the relative vorticity may certainly approach or exceed 
the value of the Coriolis parameter. In the oceans there are innumerable 
submarine features, covering a spectrum of size scales, that rise well above 
the mean surrounding depths. These observations on the inherent limitations 
of the perturbations techniques, coupled with the well known importance of 
the variation of the Coriolis parameter with latitude, provide the basis 
for the present investigation. 



8 



In subsequent chapters, Gbrtler's (19U1) conclusions will be 
accepted and it will be assumed that the flow of the barotropic ocean is 
frictionless, although it is recognized that friction must be important 
in some limiting sense. Indeed, Gttrtler pointed out that frictional effects 
may become important when the horizontal scale of the topographic feature 
under consideration is very much larger than the inertial radius of the 
flow. As the objective of this investigation is to study the combined 
effects of the variation of Coriolis parameter and finite changes in depth 
over a range of scales, subsequent assumptions will not restrict these 
essential features of the problem. 

The concept of vcrticity has become fundamental to the investigation 
of both atmospheric and oceanic flows. In its general form, the vorticity 
equation is a concise statement of one of the laws of fluid motion. For 
the class of flow under consideration its form facilitates investigation 
from either an Eulerian or Lagrangian viewpoint. This property of the 
vorticity equation will be utilized in succeeding chapters. 

The deep currents in the oceans will be treated as though they were 
barotropic, not because this mode is necessarily strictly representative of 
in situ deep currents, but because the barotropic equations of motion are 
simpler to treat exactly than are the baroclinic equations. The influence 
on this barotropic flow of topographic features with the scales of features 
existing in the real ocean will be investigated. 

Topographic effects observable in tne baroclinic flow (see Figure 1) 
must be related to those occurring in the barotropic flow, but this relation 
has not yet been determined. 
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CHAPTER II 

FORMULATION OF THE PROBLEM 



II-l. 3aslc Equations 

We consider the horizontal components of the equation of motion of 
a frictionless, homogeneous ocean or atmosphere. 



Du 



- fv 



_ 1 <)p 



Dt " f <)x ’ Dt 
together with the continuity equation 

^u . r)v ))w 

n. ^ ’ 



— + fu ° - ~ — 



_ 1 

? dy ’ 



( 1 ) 



( 2 ) 



where u, v, and w are the velocity components along the x (eastward), 
y (poleward), and z (upward) axes. The Coriolis parameter (equal to twice 
the local vertical component of the earth's rotation) is denoted by f, and 
P is the local pressure. The local acceleration of gravity is g and (° is 
the fluid density, both of which are assumed constant. 

II-2. The Vorticity Equation 

The two components of equation (1) are cross-differentiated with 
respect to y and x respectively, and then subtracted to give, after 
rearrangement of terms, 



SLh t ) . . ( £. f) b 

Dt KI J \*x 



k>u + 3v 


+ 


Jw <)u _ ^v 


[^x 3y 




Dy <3z i>x Tz 



(3) 



where ^ ~ , the vertical component of relative vorticity. 

Equation (3) may be considered as the primitive vorticity equation, 
where the change in the absolute vertical vorticity of a parcel of fluid is 
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seen to be due to the local two-dimensional velocity divergence and the 
conversion of horizontal vorticity to vertical vorticity. 

If we now consider a and v to be independent of z, the vortex tube 
term in equation (3) drops out, and we may integrate equation (2) from the 
bottom of a column to the free surface to obtain 



where is the elevation of the free surface above its undisturbed position, 
and h is the height of the undisturbed free surface above the rigid bottom, 
or ground. Substitution of the above equation into equation (3) then yields 



the equation commonly interpreted as expressing the conservation of potential 
vorticity along a trajectory. It is of interest to note that no assumption 
of steady state has been made in the preceding development. 

II-3. The Transport Streamf unction 

To facilitate the introduction of the volume transport streamfunction, 
the vertically integrated continuity equation may be written as 



For stationary conditions this equation is satisfied identically by a volume 
transport streamfunction defined such that 



jv 1 D(h + >1 ) 

55 “ (h + ?\) Dt 




(U) 





* y(x,y). (5) 



Henceforth, we will neglect the elevation of the free surface compared to 
the depth h, both as a readily justifiable approximation and for expediency. 
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In terms of the streamfunction y , the first integral of the 



vorticity equation, equation (U), is then 

1 

7-(aVf) » f 

h 



constant (along a streamline). 



( 6 ) 



Charney (1955) and Fofonoff (1962) have shown that there must be a 
functional dependence of the Coriolis parameter and the depth h on the 
streamfunction f . This may be readily seen by considering that the above 
equation is applicable along an arbitrary streamline, and that the choice 
of a specific streamline then determines the geographic path followed by a 
fluid column and hence its future depth and latitude. Substantial 
analytical difficulties arise when an attempt is made to express both the 
bathymetric relief and the variable Coriolis parameter as functions of the 
upstream flow pattern. 
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CHAPTER III 

SOLUTION FOR ZONAL FLOWS 

The difficulties mentioned at the end of Chapter II may be alleviated 
by restricting our consideration to flows that are strictly zonal upstream 
of the region of interest, and to topography that is independent of y. 

III-l. Simplification of the Problem 

Specifically, we now establish the following conditions: at some 

point upstream of the region of interest, say at x * Xq, the flow is 
defined by Y * - UHy, where H is a constant depth and U is the transport 
velocity (positive eastward )j in the first region of interest, say x^Xq, 
the depth is a function of x only, h(x). In all regions the Coriolis para- 
meter is proportional to y, the coefficient of proportionality being /?, and 
with the origin of the y axis at the equator, f * /?y. Denoting the upstream 
region by the subscript 1, we then have f j * -UHy}, hence fi(^}) ■ 

/5 yi ( f i) a - /9fj/UH. Substituting into equation (6), 

_ . y. (5 W) * ly ^ (?) 

UH 2 h W 

where V ■ (because equation (6) applies along a streamline). 

Rearranging and dropping the now superfluous subscript 1, 



(6) 



h(x)'^7 2 ^ - jy + y = - h 2 (x)/ffy . 

dx <)x UH 2 



i 
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III-2. Matching Conditions 

The flow defined by equation (8) must satisfy continuity requirements, 
that is, it must match the upstream transport at x ■ x Q . As the upstream 
transport is geostrophic, we might meet these requirements by matching 
meridional pressure gradients, but this is found equivalent to directly 
matching total transport, or transport components. The matching conditions 
are then 



- o 

c -)X dx 



dl = » - UH 

J y e) y 



at x 



(9) 



III-3 . The General Solution 

The vorticity equation (3) may readily be reduced to an ordinary 
differential equation by utilizing the separation ^(x,y) * y 0(x) . 

Upon substitution and cancellation of y, equation (8) becomes 



h(x)0" (x) - h'(x)/'(x) 




0 



h d S 



( 10 ) 



where p 2 = . 

The matching conditions (9) become 



0U O ) 

0(>-o) 



- UH 
0 



( 11 ) 



To facilitate the solution of (10) subject to (11) we make the 



following transformation of the independent variable x: 

rx 






Let 0(x) ** G(s), where s * nj h(x)dx, and s Q ■ s(x 0 ) ■ 0, 



and define D(s) by: h(x) - - UH 2 D(s). 
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Upon substitution, equation (10) becomes 

G 2 + G 2 - l/D 2 (s) (12) 

and tne matching conditions (12) becomes 

G 2 (0) - G 1 (0) - - UH 

, , U3) 

G 2 (0) - G 1 (0) - 0 , 

where the subscript 2 indicates the second region considered, that is, the 
region immediately downstream of the uniform zonal flow. Anticipating that 
we may desire to consider more than one downstream region, we may 
immediately generalize equations (12) and (13) to: 

G„ + G n * 1/D n (s) (1 h) 



and 



G n< s n „> 





2, 3, h, ....) (15) 



where 



s(a) 




h(x)dx 



The complete solution of equation ( lU ) is 



'NSr 



G n (s n ) = A n cos s n + B n sin s n + / sin(s n - t) dt 



(16) 



3 no 



D n (t) 



wnere the coefficients A n and B n are to be determined from equations (13) 
and ( 15 ) • 

The solution (16) may be evaluated (if only numerically) for any 
bathymetric profile that is independent of y. It represents the topographic 
disturbance of a zonal current subject to the following assumptions: 



1 . 
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The ocean (or atmosphere) may be idealized as homogeneous with 
a free upper surface. 

2 . Horizontal velocities are invariant with depth. 

3. The influence of friction is negligible. 

It. The convergence of the meridians may be neglected, and the 
meridional variation of the Coriolis parameter may be 
approximated as a linear function of distance from the equator 
(5-plane approximation). 

h. Displacements of the free surface from a level surface are small 
compared to the total height of the fluid column. 

6. The initial flow is zonal with zero relative vorticity. 

7. The flow, and hence the pressure field, is stationary. 



